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Motivated by the accumulating hints of new sterile neutrino species at the eV scale, we explore 
the consequences of such an hypothesis on the solar sector phenomenology. After introducing the 
theoretical formalism needed to describe the MSW conversion of solar neutrinos in the presence 
of one (or more) sterile neutrino state(s) located "far" from the (v\, v-z) "doublet", we perform a 
quantitative analysis of the available experimental results, focusing on the electron neutrino mixing. 
We find that the present data posses a sensitivity to the amplitude of the lepton mixing matrix 
element U e 4 — encoding the admixture of the electron neutrino with a new mass eigenstate — 
which is comparable to that achieved on the standard matrix element f/ e 3- In addition, and more 
importantly, our analysis evidences that, in a 4-flavor framework, the current preference for |f/ e 3| 7^ 
, is indistinguishable from that for \U^\ 7^ 0, having both a similar statistical significance (which is 

■ ~ 1.3(7 adopting the old reactor fluxes determinations, and ~ 1.8<r using their new estimates.) 

We also point out that, differently from the standard 3-flavor case, in a 3+1 scheme the Dirac 
CP-violating phases cannot be eliminated from the description of solar neutrino conversions. 
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PACS numbers: 14.60. Pq, 14.60.St 



I. INTRODUCTION 



(-} I The interest around the possible manifestation of novel neutrino properties in very short baseline (VSBL) setups has 
. been recently reawakened by the emergence of new (and the reappraisal of old) inconsistencies. Improved calculations 
' of the reactor antineutrino spectra [l[ suggest fluxes which are ~ 3% higher than previous estimates 0-Q , and have 
raised the so-called antineutrino anomaly [6|, consisting in a deficit of (almost) all the VSBL reactor measurements 
' performed at distances L < 100 m. An apparently unrelated deficit had been already evidenced in the calibration 
measurements conducted at the solar neutrino experiments GALLEX and SAGE 0, [1[ , employing radioactive sources 
' placed inside the detectors. Such two discrepancies join those already recorded at the VSBL accelerator experiments 0- 
' llll j adding further confusion to an already intricate "VSBL neutrino puzzle" . 

While some of the (old and new) anomalies may be imputable to theoretical and/or experimental inaccuracies, 
the possibility exists that they may represent a manifestation of new physics. From this perspective, they seem to 
, point towards a phenomenon of (dis-)appearance of the electron neutrinos, possibly mediated by oscillations into 
new sterile specie(s) @, Such an hypothesis finds support in the lesser degree of tension — although still a 
considerable one — now existing between the (weakened) VSBL reactor limits and the (unchanged) VSBL accelerator 
constraints, within all schemes endowed with additional sterile sjjecies [HI, [l3jj. Notably, an independent hint in the 
same direction arises from the latest cosmological data analyses [14[ , which favor a sig nificant extra relativistic energy 
content, although pointing towards masses of the new presumptive light particles |l4L[l5j. which are at the borderline 
of those (larger) suggested for the sterile neutrinos by the oscillation data. 1 Such a trend is partially corroborated by 
the latest constraints coming from primordial nucleosynthesis (Tel [l9j , which can easily accommodate one, but hardly 
two additional sterile species. 

The putative sterile neutrinos must be introduced without spoiling the basic success of the standard 3-flavor 
paradigm. This can be achieved in the so-called 3 + s schemes, where the s new mass eigenstates are assumed 
to be separated from the three standard ones by large splittings, 2 giving rise to the hierarchal pattern 3 |m| — 
m\\ <C |m§ — m\\ <C \ra 1 k — m\\ (k = 4, ...,3 + s), which ensures that the fast oscillations induced by the new mass 



1 It must be stressed that the degeneracy between the neutrino masses and the dark energy equation of state parameter w Il6| may 
allow the accommodation of such bigger masses at the cost of assuming w < — 1 and a lower expansion age of the universe, as recently 
remarked in |17| . 

2 Another (possibly coexisting) realization is provided by new sterile species separated from the active ones by extremely small mass- 
squared splittings [20I. |21||. However, such schemes do not have observable effects in VSBL settings and we do not consider them in this 
work. 

3 The solar data are insensitive to the reciprocal ordering of the third and fourth mass eigenstates, as well as to their ordering respect 
to the "doublet" (1/1, ^2). In this work we assume for definiteness mi < 7712 < m^ < 7714, always referring to the positive mass-squared 
splittings defined as Am^ = m 2 — m 2 (i < j). 
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cigcnstates are completely averaged in all settings sensitive to the Amf 2 -driven (solar) and Am^-driven (atmospheric) 
transitions, leaving unmodified the two well-established oscillation frequencies. With the additional assumption of a 
small admixture of the active flavors with the new mass eigenstates, the 3 + s schemes leave basically unaltered also 
the standard oscillation amplitudes, thus realizing a genuine perturbation of the leading 3-flavor scenario, whose size 
must respect the constraints imposed by all the existing phenomenology. 4 

The solar sector data (Solar and KamLAND) have had a pivotal role in establishing and shaping the 3-flavor 
framework, and continue to be of extreme importance in sharpening its basic parameters. As a matter of fact, 
these data are the only ones sensitive to the admixture of the electron neutrino with the first two mass eigenstates 
(vi, 1/2), also possessing a subleading sensitivity to the (averaged) 1/3-driven oscillations. This translates into stringent 
constraints on the amplitude of the elements (U e x, U e 2, U e ^) entering the first row of the lepton mixing matrix. 
Therefore, it is of certain interest to explore how such a sector "responds" to the perturbations generated by a 
non-zero mixing of the electron neutrino with new sterile species, as hinted at by the recent VSBL findings. 

The impact of new sterile species on t he p henomenology of the solar neutrino sector has been investigated in several 
works, both as a leading mechanism [2tll27| (in the "pre-KamLAND era") and as a subdominant one (2814341 ] (after the 
KamLAND results). However, all the existing analyses 5 have been performed in the simplified framework — whose 
formalism was originally developed in [35j — of pure (y\-V2) -driven oscillations, which neglects the possible mixing of 
the electron neutrino with the third standard mass eigenstate (C/ e 3 = 0) and with a new fourth one (t/ e 4 = 0). In the 
past, both assumptions were justified by the limited sensitivity of this dataset, and by the strong upper bounds put 
on these matrix elements by the reactor experiments performed with short |36j (sensitive to J7 e 3) and very short |37j 
(sensitive to U e ±) baselines. Moreover, no hint of transitions into sterile states was evidenced in the aforementioned 
analyses, so there was no reason to extend them beyond such a simple scheme. However, this situation has gradually 
changed in the recent years. On the one hand, we have witnessed a substantial increase in the sensitivity of the solar 
sector to possible departures from the simple 2-flavor approximation; the recent indication of #13 > [38j and the 
possible hint of non-standard MSW dynamics [H| (see also [40|) suggested by these data testify such a new trend. 
On the other hand, the new anomalies directly point towards a relatively big amplitude of £/ e 4, which should be now 
testable at the solar sector. 

The incorporation of the 3 + s schemes in the description of solar neutrino conversions is not a trivial task, however, 
as it requires the treatment of the MSW effect in the presence of sterile species. This problem has been recently 
addressed in [4l| . where a parameterization independent form of the lepton mixing matrix has been exploited, and 
selected numerical examples have been given for the relevant transition probabilities. In this work we make a step 
forward and, by adopting a convenient parameterization of the mixing matrix, we put quantitative constraints on 
these schemes, focusing on the electron neutrino mixing. The rest of our paper is organized as follows. In Sec. II 
we introduce the 3+1 neutrino framework and present the basic formulae needed to interpret the flavor oscillations 
in vacuum and in (solar) matter within such a scheme. In Sec. Ill we discuss the results of the numerical analysis, 
drawing our conclusions in Sec. IV. Four appendices address the following (more technical) issues: A) The treatment 
of the solar MSW effect in a 3 + 1 scheme; B) Its generalization to the 3 + s frameworks; C) The incorporation of 
Earth-induced matter effects; D) The potential sensitivity of solar neutrino flavor transitions to the Dirac CP-violating 
phases entering the 4-lepton mixing matrix. 

II. FRAMEWORK AND BASIC ANALYTICAL RESULTS 
A. Parameterization of the mixing matrix 

In the presence of a fourth sterile neutrino v S7 the flavor (y a , a = e, n, r, s) and the mass eigenstates (i/j, i = 1,2, 3, 4), 
are connected through a 4 x 4 unitary mixing matrix U, which depends on six complex parameters [42j. Such a matrix 
can thus be expressed as the product of six complex elementary rotations, which define six real mixing angles and six 
CP-violating phases. Of the six phases three are of the Majorana type and are unobservable in oscillation processes, 
while the three remaining ones are of the Dirac type. For simplicity, in this work, we set to zero all the Dirac phases, 
commenting only in Appendix iDl on the potential sensitivity of the solar data to them. 

As it will appear clear in what follows, for the treatment of the solar MSW transitions under study, it is convenient 



4 This includes the atmospheric neutrino data collected at neutrino telescopes [23l . where the new sterile species may leave distinctive 
imprints [HQ. 

5 An exception is constituted by the work [30| . where a qualitative discussion of the impact of non-zero U e 3 is provided. 
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to parameterize the mixing matrix as 

U = R23R24R34R14R13R12 = BR14R13R12 = AR\2 j (1) 
where Rij represents a real 4x4 rotation in the plane containing the 2x2 submatrix 



Rt* z = ^ 13 , (2) 

\ s ij ''■ i J 

in the subblock, with (dj = cos 8ij,Sij = sindij), and the matrices A = UR\ 2 and B = R23R24R34 have been 
introduced for later convenience. 

The parameterization in Eq. ([1]) has the following properties: I) For vanishing mixing involving the fourth state 
(614 = 824 = $34 = 0) Eq. (p} reduces to the 3-flavor mixing matrix in its standard parameterization [4j|; II) The 
leftmost positioning of the matrix i? 2 3 allows us to rotate away the mixing angle #23 from the solar MSW dynamics 
(see Appendix [5J and from the expression of the mixing elements involving the sterile flavor (see the discussion below); 

III) The positioning of the product R14R13 close to the rightmost matrix i?i 2 makes the corresponding mixing angles 
appear in a symmetrical way in the expressions of the admixtures of the electron neutrino with the (y 1,1/2) mass 
eigenstates, inducing small admixtures with the "far" (^3, V4) mass eigenstates, which, in the limit of small #14 also 
appear in a symmetrical form, being £/ 2 3 ~ sin 2 613 and £/ 2 4 = sin 2 614. In fact, the mixing matrix elements involving 
the electron neutrino are expressed as 

U e i = An C12 - A 12 S12 , (3) 

Ue2 = An S\2 + A12 C12 , (4) 
Ue3 = A 13 , (5) 
Uei = A u , (6) 

in terms of the elements of the first row of the matrix A. The rightmost positioning of the product R14R13R12 in such 
a matrix ensures that its element A12 is equal to zero, thus leading to the explicit expressions 

(7) 
(8) 
(9) 
(10) 

It should be noted that these expressions are not affected by the pre- multiplication of the matrix B in Eq. ([IJ , since 
this induces a rotation in the hyperplane orthogonal to the axis with index i = 1, and therefore are independent of the 
specific order of the three matrices (-R23, -R24, -^34)- Since the solar data are sensitive also to transitions into sterile 
neutrinos, we will need the expressions of the matrix elements involving the sterile flavor, which, taking into account 
the last equality in Eq. ([1]), can be written as 

Usi = A41 C12 - A 42 S12 , (11) 
U s2 = A 41 s 12 + A 42 C12 , (12) 
U s3 = A i3 , (13) 
U sA = A 44 , (14) 

in terms of the elements of the fourth row of the matrix A, which, in our parameterization have the explicit expressions 

Ail = C 2 4(S34S13 - C34S14C13) , (15) 

A}2 = -S24 , (16) 

-443 = -C24(S34C13 + C34S14S13) , (17) 

-444 = C24C34C14 , (18) 

independent of the mixing angle #23, due to the leftmost positioning of the matrix i?23 in the definition of the mixing 
matrix in Eq. ((T|). The relations in Eq. (jlllllSp also evidence the following further features of our parameterization: 

IV) For small values of all the mixing angles involving the fourth mass eigenstate [s 2 4 < l(i = 1,2,3)], the sterile 
flavor content is mostly distributed on the fourth mass eigenstate (U S 4 ~ 1); V) For # 13 ~ 6* 14 ~ 0, the sterile content 
of the first two mass eigenstates essentially depends only on (9 2 4, being in this limit + t/ 2 2 — -4 2 2 = s 2 4 ; VI) For 
813 ~ 814 ~ 824 — 0, the mixing angle #34 basically exchanges the amplitude of U S 3 with that of U S 4. 



U e l 


= C14C13C12 


Uei 


= C14C13S12 


Ue3 


= C14S13 , 


Ue4 


= S14 . 
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It should be stressed that our parameterization — "tailored" for the solar sector — is different from that commonly 
adopted for the 4- flavor analyses of the atmospheric and long baseline neutrino data [3|4^]. In such a case, the mixing 
matrix is still taken of the form U = BR14R13R12, but with a different choice of the ordering of the rotations entering 
B (B atm = R34R24R23) 1 which is dictated by the fact that — > v T transitions are sensitive to the matrix elements 
connecting and v T to the new mass eigenstate z/4. Indeed, with such a choice, in the limit of small admixtures 
with the fourth state, one has the approximate expressions s\ A ~ and s§ 4 ~ L^ 4 . In this regard, we observe that 
the parameterization we have adopted could be classified as of a "mixed form" , in between the two ones recognized 
in |33[ as more apt for studying those datasets possessing, respectively, a prevailing sensitivity to admixtures (of the 
new state v s ~ 1/4) with the flavor eigenstates {U a A 7^ 0) or with the mass ones (U s i 7^ 0). Such a particular form 
ensues from the peculiar properties of the solar sector data, which are sensitive both to U e 4 ( "favor- type" admixture) 
and to the U s i's ("mass- type" admixtures). 



B. Four-flavor evolution 

The evolution of the neutrino flavor eigenstates is governed by the Schrddinger-like equation 



.d_ 

dx 



( v e 



2 1=^/1 £1. < 19 ) 



where the Hamiltonian 



H f = Hf n + Uf n = UKU T + V(x) , (20) 



has been split in the sum of a kinematical and a dynamical term. In Eq. (|20| K denotes the diagonal matrix containing 
the wavenumbers ki = mf/2E(i = 1,2,3,4) (rrii and E being the neutrino mass-squared and energy respectively), 
while the matrix V(x) incorporates the matter MSW potential 0, |4tJ- Barring irrelevant factors proportional to the 
identity, we can define the diagonal matrix containing the three relevant wavenumbers as 

K = diag(0, k so i,k atm , k new ) (21) 

/ Ami Am? n Am?, \ , . 

and the matrix encoding the matter effects, as 

V = diag(y C c, 0, 0, -V NC ) , (23) 

where 

Vcc = V2G F N e (x), (24) 

is the charged-current interaction potential of the electron neutrinos with the background electrons having number 
density N e , and 

V NC = ~V2G F N n {x), (25) 

is the neutral-current interaction potential (common to all the active neutrino species) with the background neutrons 
having number density N n . For later convenience, we also introduce the position-dependent parameter r x = r(x) 
defined as the positive-definite ratio 

_ V N c(x) _ 1 N n (x) 

V CC (x) 2N e (x)- ( ' 
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C. Oscillation probabilities in vacuum 

In the case of propagation in vacuum, Eq. (|19[) leads to the survival probability of electron (anti-)neutrinos (which 
is relevant for the reactor experiment KamLAND 6 ) 

P ee = 1-±Y. U ej U ek sin 2 , (27) 

j>k 

with the elements U e i (i = 1, 4) determining the amplitudes of the oscillating terms having phases 4>jk — Arrij k L / %E 
developed over the baseline L. In the hierarchical limit k so i <C k atm <C k new the fast oscillations induced by the two 
larger wavenumbers are completely averaged and the survival probability can be written as 

Pee = (1 " E U ekf P e V e + £ Utk > (28) 
fc=3 fc=3 

where 

P^ = l-4s 2 2 c? 2 sin 2 <£ 12 , (29) 

is the well known 2-flavor expression of the survival probability in vacuum. Eq. (|28|) shows that the presence of the 
"far" eigenstates v% and is felt as a lack of unitarity of the (1*1,1*2) sector and that there is an exact degeneracy 
between C/ e 3 and U e i- In our parameterization of the mixing matrix Eq. (|28[) reads 

Pee = c 14 C 13^ee + C 14 S 13 + S 14 > (30) 

which implies an approximate degeneracy between the two (small) mixing angles #13 and #14. 

D. Transition probabilities in matter 

Matter effects play a central role in the conversion of solar neutrinos and should be incorporated following the 
treatment presented in Appendix |A"1 It must be observed that, although the solar data are mainly sensitive to the 
survival probability of the electron neutrinos, they also posses a sensitivity to the transition probability into sterile 
states (Pes) through the neutral current (NC) measurements performed by the SNO experiment and, to a lesser 
extent, through the elastic scattering (ES) interactions exploited by SuperKamiokande and Borcxino. 

As shown in Appendix [A] for the small values of the two mixing angles #13 and 9u we are considering, the 
propagation of solar neutrinos is adiabatic, and the transition probabilities only depend upon their production and 
detection points. Neglecting Earth- induced matter effects 7 one has the general expressions 

4 

%^^)-^^(C) 2 (a = e,fi,T,s), (31) 

i=l 

where the [/™'s denote the mixing matrix elements involving the electron neutrino in the production point. As shown 
in Appendix A, in our parameterization these elements can be simply obtained by replacing in Eqs. (|71 ll0[) the mixing 
angle #12 in vacuum with the corresponding one in matter 



U n e l = c 14 c 13 c? 2 , (32) 

U"i = C14C1 3 S" 2 , (33) 

US = Ue3 = C14S13 , (34) 

Uel = UeA = SiA , (35) 



In this discussion we neglect the small matter effects induced by the interaction of the electron antineutrinos with the Earth crust. 
However, for the sake of precision, we include these effects in the numerical analysis presented in the next section. 

Here, for simplicity, we neglect Earth matter effects which, however, are properly included in the numerical analysis following the 
prescription described in Appendix [Cl 
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where 0™ 2 can be calculated using the prescriptions provided in Eqs. (|A18tiA20|) . In the numerical analysis we are 
going to present in this work, we will limit ourselves to the simple case #24 = $34 = 0. In this case, the matrix A in 
Eq. ([I]) takes the simpler form A = R23R14R13, implying the following simplified expressions for the mixing elements 
involving the sterile flavor 



Usl = -S14C13C12 , 

U S 2 = -S14C13S12 , 

U S 3 = -S14S13 , 
U S 4 = C14 , 



(36) 
(37) 
(38) 
(39) 



which show that, in this case, the small departure from unity of the sterile flavor content of the fourth mass eigenstatc 
is redistributed (by unitarity) 
one arrives at the expressions 



is redistributed (by unitarity) essentially to the first two mass cigenstates, leaving U 2 3 ~ 0. Using Eqs. (|A18IIA20|) . 



'14 c 13^ee + c 14 s 13 



'14 ' 



P - r 2 s 2 r 4 P 2v 
r es — c 14*14 L 13- r ee 



2 2 4 

c 14 s 14 s 13 



where 



2 / m\2 



(40) 
(41) 

(42) 



is the well known expression of the 2-flavor survival probability in the adiabatic regime [48| but with the mixing angle 



6*12 obtained using the rescaled charged-current potential Vcc ~ > ( c 
of #14 = 0, Eqs. (|40H4ip return the standard 3-flavor expressions [49| 



p3u _ 4 
r ee ^ c 13 

P Zv = 

± es u ' 



P. 



1v 



'13 ! 



with the well known rescaling of the standard MSW potential Vcc —> c\ z Vcc [49LI5C 
0i 4 7^ 0), one gets 



2c 



Pee (013 - 0) = C 4 4 ^ 

P es (0 13 = 0) = c 2 u sUP 2 » + 1) 



■;Si 4 )c{ 3 Vcc (see Appendix A). In the case 

(43) 
(44) 

Conversely, in the case (#13 = 0, 

(45) 
(46) 



with the rescaled potential Vcc (c 2 4 + r x s 2 4 )Vcc ■ We observe that, while the form of the electron neutrino survival 
probability is identical to that obtained in the 3-flavor case, modulo the replacement #i 3 — > On and a small change in 
the rescaling factor of the potential, in this case, there is a non-zero transition probability into sterile states. In all 
cases, the change in the potential is very small and, analogously to the 3-flavor limit (5l| . it introduces only a mild 
(undetectable) energy dependence of the (significant) energy-independent suppression induced by the kinematical 



facto 



r c 14 c 13 



1 



ZA 14 



2sf 3 appearing in Eq. fl40 



III. NUMERICAL RESULTS 

In our analysis we have included the data from the Chlorine (CI) experiment [52|, the Gallium (Ga) detectors 
SAGE HI and GALLEX/GNO (H-H^, Super-Kamiokande H3 (SK), all the three phases of the Sudbury Neutrino 
Observatory [H-[6l| (SNO), and Borexin o l62l l63j (BX). We have not included the spectral information provided by 
the SNO Low Energy Threshold Analysis [64| , since this can be used only under the assumption of unitary conversion 
among the active neutrino species (P ee + Pen + Per = 1), which is violated in the presence of transitions into new 
sterile states. For the sake of precision, we have also incorporated the small Earth-induced matter effects, following 
the treatment presented in Appendix [C] Concerning KamLAND, we have included in our analysis the latest data 
released in (65| . For definiteness, we have adopted the new improved reactor flux determinations All plots will 
refer to this case and, when appropriate, we will comment in the text on the differences that would arise with a 
different choice of the reactor fluxes. In all numerical computations we have set 6*24 = #34 = 0, commenting at the 
end of the Section on the possible role of these parameters. Therefore, the parameter space spanned by our analysis 
will involve the solar mass-spliting Am 2 2 and the three mixing angles (#12, #13, On)- 

We start our numerical study considering the familiar three-flavor case (6*i 3 7^ 0,#i4 = 0), in which the results of 
the analysis depend on the three parameters (Am 2 2 , #12, #13)- This case will serve as a useful term of comparison for 
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FIG. 1: Region allowed in the [sin 2 #12, sin 2 613] plane for #14 = 0, after marginalization of Am 2 2 as constrained by KamLAND, 
separately (left panel) by solar (S) and KamLAND (K) data and by their combination (right panel). In both panels it has been 
set 824 = #34 = 0. The contours refer to A\ 2 = 1 (dotted line) and A^ 2 = 4 (solid line). 



the more general results of the Av analysis. In the left panel of Fig. 1 we show the region allowed by solar (S) and 
KamLAND (K) in the plane spanned by the two mixing angles, having marginalized away the solar mass splitting 
in the region determined by KamLAND. Respect to previous analyses [38l. I65l469l | . the KamLAND data taken alone 
now tend to prefer values of #13 > (see also 0, [7(|). This behavior can be traced to our adoption of the new 
(higher) reactor fluxes [l|. In fact, according to Eq. (pIB"]) . a larger value of #13 is now required to sup press the bigger 
total rate induced by the new higher fluxes. Furthermore, similarly to previous analyses (38l. [6Fj| - l69j . for #13 > the 
values of the mixing angle 612 identified by the solar and KamLAND experiments are in better agreement due to the 
opposite-leaning correlations exhibited by their respective contours, giving rise to an enhanced preference for non-zero 
6*i3 in their combination (right panel). We find that the 2-flavor case (6*13 = 0) is disfavored at the 1.8 a level (which 
is reduced to 1.3 a using the old reactor fluxes). 

As as second step we switch on only the mixing angle #14, setting #13 = 0. In this case, the results of the analysis 
depend on the three parameters (Am^ 2 , #12, #14), whose allowed regions are displayed in Fig. 2. As discussed in 
Sec. II, KamLAND cannot distinguish #13 from #14 and, as a result, the region identified by such an experiment is 
identical to that found in the 3-flavor case. In contrast, the region determined by the solar data is slightly different 
from the corresponding one identified in the 3-flavor case. In particular, we see that the correlation in the [s^Su] 
plane is different from that exhibited in the [s\ 2 , S13] plane. To understand this point, it is useful to write the relations 
connecting the total fluxes observed in the solar neutrino experiments to the electron neutrino survival probability 
and their transition probability into sterile states. Taking into account that P ee + P efi + P eT + P es = 1, we have 

4> cc = $ B (P ee ), (47) 
ES = $ s ((P ee ) + rv(l - (P ee ) - (P es )) , (48) 
NC = <Ml-(P es )), (49) 

where $b is the solar 8 B neutrino flux expected in the absence of oscillations, while </> x (X = CC, ES, NC) is the flux 
expected in charged-current (CC), elastic-scatte ring (E S) and neutral-current (NC) reactions. The symbol () denotes 
the average over appropriate response functions |7ll . l72| , which depend on the experiment and on the specific reaction 
considered. We remind that: I) The radiochemical experiments (CI, Ga) employ a CC absorption reaction and are 
thus sensitive only to P ee ; II) The BX and SK detectors make use of an ES reaction and can probe also P es , although 
with sensitivity reduced by the small ratio r a ~ 1/6 of the energy-averaged ES cross-sections of v^ T and v e ; III) The 
SNO experiment detects the solar neutrinos using all the three reactions (CC, ES, NC), thus possessing a pronounced 
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FIG. 2: Region allowed in the [sin 2 #12, sin 2 On] plane for #13 = 0, after marginalization of Am 2 2 as constrained by KamLAND, 
separately (left panel) by solar (S) and KamLAND (K) data and by their combination (right panel). In both panels it has been 
set 824 = 834 = 0. The contours refer to A\ 2 = 1 (dotted line) and A\ 2 = 4 (solid line). 



sensitivity to P es through the NC process. 

In the 3-flavor case, the sensitivity of the solar data taken alone to the mixing angle #13 arises from an interplay of 
low-energy (LE) and high-energy (HE) (essentially the SNO CC /NC ratio) data [73, [74[ , engendered by a different 
dependence of the survival probability on the two mixing angles #12 and #13 in the two regimes. From Eqs. (|47ti49|) . 
taking into account that in the LE (vacuum-like) regime ~ 1 — %s\ 2 c\ 2 and in the HE (matter-dominated) regime 
~ s\ 2 , and using Eq. (j4~5T) . we have the following relations for the LE and HE fluxes 

K pLE ^ (1 _ 2s?a)(1 _ 2s 2 2C 2 2) m 

~ P™^(1-2 S ? 3 ) S ? 2 . (51) 

HE 

The different relative sign of the two factors proportional to sf 3 and s\ 2 m Eqs. (|50ti51|) gives rise to an opposite 
correlation among these two parameters, providing an enhanced sensitivity to small departures from zero of the 
mixing angle #13. Since the SNO CC/NC ratio is measured with better precision than the LE flux (essentially 
provided by the Ga experiments), the negative relative sign in Eq. (|51|) prevails in the global fit, giving rise to the 
positive overall correlation among the two mixing angles appearing in the left panel of Fig. 1 in the curve designed 
with label "S". Indeed, according to Eq. (|5Tj) , a bigger value of #13 is needed to counterbalance the effect of a larger 
6*12, in order to keep the CC/NC ratio at the fixed value determined by the SNO experiment. 

In the case of On ^ 0, while we have an identical expression for the expected flux in the LE limit (modulo the 
replacement #13 — > On), the SNO CC/NC ratio depends also on P es . For small values of On, using Eqs (|45H46[) . we 
have 

^ ex P e LE ~ (1 - 2 S 2 4 )(1 _ 2s \ 2 c\ 2 ) , (52) 
c P e « E (l+P e ™)^(l-c? 2S ? 4 ) S ? 2 . (53) 

HE 

In this case the coefficient in front of the product s^ 4 s 2 2 in the SNO CC/NC ratio is still negative but three times 
smaller (c\ 2 ~ 0.7) than the corresponding one appearing in the 3-flavor expression. As a result, in such a case, the 
synergetic effect of the combination LE and HE data slightly decreases, with a lower sensitivity to On and a weak 
negative overall correlation in the [sl 2 ,sl 4 ] plane (see the curve designed with label "S" in the left panel of Fig. 2), 
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Solar + KamLAND 




0.02 0.04 0.06 0.08 0.1 



sin 2/ $ 13 



FIG. 3: Region allowed, after marginalization of Am 12 and #12, by the combination of solar and KamLAND data (for 024 = 
O34 = 0). The contours refer to A\ 2 = 1 (dotted line) and A\ 2 = 4 (solid line). 

as now the LE data can compete with the HE ones in determining the relative sign in front of the product s^ 4 sf 2 - 
Similarly to the 3-flavor case, the values of the mixing angle 812 identified, respectively, by solar and KamLAND are in 
better agreement for #44 > 0, with an enhanced preference for non-zero values of this parameter in their combination 
(right panel in Fig. 2). Also in this case we find that #44 = is disfavored at the 1.8 a level (which is reduced to 1.3 a 
using the old reactor fluxes). 

The following small differences appear between the two cases: I) A weaker upper bound on #14 (sf 4 < 0.089 at the 
2cr level) with respect to that obtained for 6*13 (sf 3 < 0.070 at the 2er level); II) A slightly bigger best fit value for 
#14 ( s i4 = 0.041) with respect to that obtained for 6*13 (s1 3 = 0.033). It is interesting to note that the best fit value 
obtained for 9u practically coincides with that indicated by the VSBL reactor and Gallium calibration anomalies 
taken in combination Q . Therefore, combining the solar sector results with such data would reinforce their preference 
for non-zero 9u, providing an overall indication, which we roughly estimate to be around the ~ 4cr level. 

As a third step of our numerical analysis, we have switched on both mixing angles (6*43 7^ 0, 6*44 7^ 0). In Fig. 3 we 
show the region allowed by the combination of solar and KamLAND in the plane spanned by such two parameters, 
having marginalized away both the mass splitting £\m\ 2 an d the mixing angle #42 . From this plot we see that there is 
a complete degeneracy among the two parameters. In practice, this dataset is basically sensitive to the combination 
Z7g 3 + C/g 4 , the small deviations from this behavior being induced by the SNO NC measurement. Therefore, the solar 
sector data, while indicating a weak preference for non-zero mixing with the "far" eigenstates v% and 1/4, cannot 
distinguish between them. The same conclusion holds true if additional sterile neutrinos are considered, as shown in 
Appendix B. 

We conclude this section with a final remark. As in our analysis we have set 624 = #34 = 0, the question arises as 
to whether a non-zero value of such parameters can alter the basic conclusions of the analysis. With this purpose, we 
have performed the analysis considering this possibility and we have found that, although the correlations observed 
in Figs. (1,2) can slightly change, the degeneracy among the two parameters $43 and #44 persists, as one may have 
expected, since the additional freedom given to the system can only worsen its preexisting degeneracies. We postpone 
to a future work the investigation of the constraints attainable on the mixing angles 624 and $34 , noticing that in this 
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case one should take into account the (small) effects of the CP violating phases entering the lepton mixing matrix, as 
discussed in Appendix D. 

IV. CONCLUSIONS 

Motivated by the recent experimental findings which point towards the existence of new light sterile neutrino species, 
we have explored their impact on the solar sector phenomenology. Working in a CPT-conserving 3+1 scheme, we 
have considered the perturbations induced by a non-negligible mixing of the electron neutrino with a fourth sterile 
neutrino specie. Our quantitative analysis shows that the present data posses a sensitivity to the amplitude of the 
lepton mixing matrix element U e ±, which is comparable to that achieved on the standard matrix element U e 3. In 
addition, our analysis evidences that, in a 4- flavor framework, the current preference for \U e 3\ is indistinguishable 
from that for \U e i\ ^ 0, having both a similar statistical significance (which is ~ 1.3<r adopting the old reactor fluxes 
determinations, and ~ 1.8a using their new estimates.) Such a degeneracy naturally extends to the more general 3 + s 
schemes, implying that in these frameworks, the present hint for non-zero #13 must be reinterpreted as a preference for 
a non-zero mixing of the electron neutrino with mass eigenstates distant from the {vi,V2) "doublet". Different kinds 
of probes are needed in order to discriminate whether such an admixture — if effectively confirmed to be non-zero — 
is realized with the third standard mass eigenstate 1/3 or with new neutrino specie(s). 
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Appendix A: The MSW effect in a 3+1 scheme 

To treat the MSW effect in a 3+1 scheme it is convenient to introduce the new basis 

v = A T v, (Al) 

where 

A = -R23-R24-R34-R14-R13 = URj 2 , (A2) 

is the same matrix defined in Eq. (TT]). In this new basis the Hamiltonian assumes the form 

H = H km + H dyn = R 12 KRj 2 + A T VA , (A3) 

and, in the hierarchical limit k so i <C k atm <C k new , similarly to the 3- flavor case (48j , one can reduce the dynamics to 
that of an effective 2v system. Indeed, from Eq. (|A3|) one has that the (3, 3) and (4, 4) entries of H , being proportional 
to katm and k new , are much bigger than all the other ones, and at the zeroth order in the small quantities Vcc /k a tm, 
VNc/kaUn, Vcc/k n ew, Vnc /knew, the third and fourth eigenvalues of H are much larger than the first two ones. As 
a result, the states P3 = 1/3 and P4 = z/4 evolve independent one of each other and, more importantly, completely 
decoupled from T>\ and v 2 . Extracting from H the submatrix with indices (1, 2) one obtains the 2x2 Hamiltonian 

H 2v = H^+H%r, (A4) 

governing the evolution of the (Pi, P2) system, whose dynamical part has the form 

fjdyn tt , s( A^ + TxA 2 ^ AuA 12 + r x A 41 A 42 \ / * r\ 

H *» = Vcc{x) Un^ + ^42 A\ 2 + r x A% ) ■ 

In our parameterization the rotation R 23 appears as the leftmost matrix in A, and therefore the product A T V A does 
not depend on #23 since -R23 commutes with the matrix V of the potential given in Eq. (|23[) . Furthermore, the matrix 



11 



A contains the product R14R13 as its rightmost factor, implying that its element Ai 2 is equal to zero. Therefore, the 
2x2 Hamiltonian in Eq. (|A5[) can be recast in the form 

frdyn T , , ( I 2 + r x a 2 r x aft \ 

H 2 * =Vcc{x)[ rxa/3 rxj32 ) , (A6) 



where the parameters (a,/3, 7) are defined as 



a = A ix = C24(S34S13 ~ C34S14C13) , (A7) 

ft = A i2 = -s 24 , (A8) 
7 = An = C13C14 , (A9) 

which, according to Eqs. (|11|14|) . are related to the mixing matrix elements as follows 

Usi = ac X2 -fts 12 , (A10) 
U a2 = asi 2 + /3ci2, (All) 

showing that the sum of a 2 and /3 2 represents the sterile content of the (ui, v 2 ) "doublet" 

a 2 + ft 2 = U 2 i + U 2 2 . (A12) 

Furthermore, using Eqs. (j5] fTU)) we have 

7 2 = U 2 i + U 2 2 = l-U 2 3 - U 2 ei . (A13) 

It is instructive to consider the following limit cases: 

I. For vanishing mixing with the fourth neutrino specie (#14 = #24 = #34 = 0) one has a = f3 = and 7 = C13, 
thus recovering the standard 3- flavor result [49|, [5(| , in which the admixture with the third state ^3 induces the 
simple rescaling of the potential Vcc —> c 2 3 Vcc, 

II. When only the mixing angle #14 is different from zero, a = — sn, (3 = and 7 = C14, with the position-dependent 
rescaling of the standard potential Vcc —> ( c i4 + v xs\a)Vcc] 



III. If both 0i3 7^ and 9n ^ (and #24 = ^34 = 0), a = — S14C13, ft — 0, 7 = C13C14, with the position-dependent 

41 



rescaling of the standard potential Vcc — > ( c i4 + r xS 2 4)c 2 3 Vc , c; 



IV. In the case of no admixture of the electron neutrino with the "far" states 1/3 and ^4 (#13 = #14 = 0), one 
has a = 0, ft = — s 2 4 and 7=1, with the position-dependent rescaling of the standard MSW potential 



Vcc (1 — r xS 2 4)Vcc , m agreement with the result found in [35| for this particular case 



In the general case the modifications are less obvious than a simple rescaling of the standard potential, as the 2x2 
Hamiltonian in Eq. (|A6[) can contain both new diagonal and off-diagonal terms. The matrix H 2l , in Eq. (|A4[) will be 
diagonalized by a 2 x 2 rotation 

*3*"m - (si 2%) • < a «> 

which defines the mixing angle in matter Q™ 2 that, in general, will depend on all the mixing angles except for #23- 
Therefore, if we define the 4x4 rotation in the (1,2) plane 

(pm 2x2 n \ 
12 /2 x 2 J , (A15) 

the starting 4-dimensional Hamiltonian in Eq. (|20[) will be diagonalized by the matrix 

U m = AR™ , (A16) 



In the different parameterization adopted in |35H for the lepton mixing matrix, the role of Sj 4 is taken by c| 3 c| 4 . In both parameteriza- 
tions, the rescaling factor has (obviously) the same physical interpretation, being determined in both cases by the sterile content of the 
(yit V2) sector. 
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which connects the flavor eigenstates to the instantaneous energy eigenstatcs in matter 



The mixing angle in matter #"2 is given by 



u r - 



(A17) 



-P sin 267k = sin 20 12 + 2v x r x a/3 , 
k 

% cos 29? 2 = cos 29 12 ~ v xl 2 ~ v x r x {a 2 - /3 2 ) , 
k 

where k = k so i and the neutrino wavenumber in matter k m is defined by 



(A18) 
(A19) 



= [cos20 12 - v xl 2 - v x r x (a 2 - f3 2 )] 2 + [sin20 12 + 2v x r x af3] 2 , 



(A20) 



with v x = Vcc{x)/k. It should be noted that for a/3 < both terms in Eq. (|A20[) may become simultaneously small, 
even for large values of #12. In particular, if at a point x along the neutrino trajectory the two conditions 9 



- = [ 1 2 + r x {a 2 ~p 2 )] 2 +ArW0 2 , 
vt 



tan 2012 = 



2r x \a{. 



7 2 +r x (a 2 - P 2 ) 



are both satisfied, the difference between the two energy eigenstates in matter approaches zero {k n 
a case one expects important non-adiabatic effects encoded by a non-zero swapping probability P c 



-+0). 



(A21) 
(A22) 
In such 



between the two energy eigenstates in matter. However, it turns out that the second condition [Eq. (|A22[) ] cannot 
be realized for realistic values of the parameters involved in the conversion of solar neutrinos. On the one hand, the 
mixing angle 6*12 is constrained to have relatively big values by the KamLAND experiment (see the left panels of 
Figs. 1-2), which sets the robust lower limit tan2#i2 > 1.0 (at the 4cr level), independent of any kind of matter-effects. 
On the other hand, the possible excursion of the right term in Eq. (|A22[) is severely constrained by the non-solar 
neutrino oscillation phenomenology and by the properties of the Sun. Indeed we have that: I) In the Sun the ratio 
r x never exceeds the maximum value (assumed at the Sun center) r max ~ 0.25; II) The mixing angles (#13, #14) are 
bounded by the reactor experiments (s 2 3 < 0.1 and s 2 4 < 0.1). These circumstances ensure that the numerator in 
the right term of Eq. (|A22|) is always small, while keeping its denominator close to one. Allowing for arbitrary values 
of the other two mixing angles (6*24,^34), we estimate that the ratio on the right side of Eq. (|A22|) never exceeds 
0.2. Therefore, non-adiabatic effects arc completely irrelevant in the problem under study, as we have explicitly 
checked by a numerical scan of the relevant parameter space. We stress that the same conclusion would not hold 
if the assumption of CPT invariance were abandoned, as it ensures that the mixing angles probed by the reactor 
antincutrinos arc identical to those involved in the conversion of the solar neutrinos. 

Considering an electron neutrino produced in the Sun, the probability to detect it on the Earth with flavor a will 

be 



(A23) 



where the U^s are the mixing matrix elements calculated in the production point, while the £, ~ Am 2 1 L/2E are the 
phases acquired by the energy eigenstates during their propagation from the Sun center to the Earth surface (L ~ 1 



Equations |)A21I IA22|| generalize analogous conditions introduced in [75| in the context of solar neutrino conversion in the presence 
of non-standard neutrino interactions of the flavor-changing type. In that study, focused on the case of small mixing angles, it was 
shown that non-adiabatic effects arc enhanced when: I) v x = cos2$i2 ~ 1 (corresponding to the common resonance condition); and II) 
tan 20i2 = — 2e e a, where the coupling constant e ea , parametrizing the strength of the new interaction between the two different flavors 
v e and v a , plays the role of our off-diagonal parameter r x afi. In the presence of additional flavor-diagonal interaction terms of the type 
e ee and too one would have the more general conditions: I) [l/v x = [(1 + t ee — e-aa) 2 + 4eg Q ]; and II) tan 2#i2 = — 2e e a/(l + fee — Coo). 
In our case the diagonal terms correspond to e ee = —1 + ■y 2 + r x a 2 and e aa = r x f} 2 [see Eq. I IA6I 1I . 
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a.u.). The information contained in the (large) phases gets lost by the spatial average over the neutrino production 
zone and by the energy smearing [76[ , and Eq (|A23[) reduces to 

4 

P(u e ^u a ) = ^2\U ai \ 2 \U^\ 2 (a = e,fi,r,s), (A24) 

which coincides with the Equation (|3"Tj) used in Sec. II. In summary, the calculation of the transition probability is 
reduced to the following steps: I) Given the four mixing angles (#13, #14, #24 7 $34), as defined in the parameterization 
in Eq. (pj, calculate the coefficients (a, ,3,7) making use of Eqs. (|A7IIA9j) : II) Determine how the mixing angle #12 in 
vacuum gets modified in matter applying Eqs. (|A18IIA20|) for the expression of Q™ 2 \ III) Deduce the electron neutrino 
mixing elements U™ in matter from Eq. (|A16|) ; IV) Derive the transition probabilities P(y e — > v a ) using Eq. (|A24I) . 



Appendix B: Generalization to a 3 + s scheme 



The generalization to more than one sterile specie is straightforward and is obtained by observing that the 2x2 
matrix H^f™ entailing the non-trivial dynamics is always given by the submatrix with indeces (1,2) of the matrix 

H dvn = A T VA, (Bl) 

where now A, still defined as A = UR[ 2 like in the 3+1 scheme, has dimension 3 + s, and the matrix V is given by 
diag{Vcc 1 0, 0, — V/vcj —Vnc, ■•■)■ If the matrix A is taken with the product Ri.3 +s ...RuRi3 of the matrices involving 
the first index as its rightmost factor (ensuring A\2 = 0), the effective 2x2 Hamiltonian H^T has the same form 
of Eq. (|A6|) provided that the (combinations of the) three coefficients (a, j3, 7) appearing in Eq. (|A6[) arc replaced as 
follows 

s 

a 2 -^^a 2 , 

i=l 

i=l 
s 

a/3 -^^ajft , 

i=l 

with 



(B2) 
(B3) 
(B4) 



at =A 3+iA (i = l,...,s), (B5) 

Pi=A w (i = l,...,s), (B6) 

while 

V-H-t^-^OS. (B7) 

i=4 

Form this last formula it is evident that the role of the electron neutrino mixing with additional sterile species is 
completely symmetrical to that played by the fourth one. 



Appendix C: Inclusion of Earth matter effects 



We briefly review the analytical results pertaining Earth matter effects, for the sake of completeness and self- 
consistency of the paper. In general, Earth matter effects intervening prior to the detection of the solar neutrinos, 
can be implemented by the following substitution in Eq. (|31l) 

Uli^Pcci (i = l,2 ! 3,4;a = e,/i,r,s), (CI) 

where P a i = P(yi — > v a ) are the conversion probabilities in the Earth of the mass eigenstates into the flavor ones. 
Although the P a i& can be calculated numerically, it is possible to simplify their evaluation by reducing the Av dynamics 
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to that of a 2v system, in analogy with the solar-matter induced effects discussed in Appendix A. In fact, the P a i's, 
can always be written as 



\{AWR 12 ) 



\(AZi 



(C2) 



where, from right to left: I) The matrix R12 rotates the initial mass eigenstates into the auxiliary flavor basis v defined 
in Eq. (|A1[) ; II) The matrix W contains the (complex) transition amplitudes among the eigenstates of such new basis, 
whose non trivial dynamics is confined to the (pi,T>2) sector; III) Finally, the matrix A rotates back the auxiliary 
basis v to the standard flavor basis. The 4x4 matrix W contains in its (1,2) subblock the non-trivial information, as 
follows 



Z 12 = WRw = 








VT - Pe""' 0\ 

VPe^ 

1 

1/ 



(C3) 



where the 2-flavor transition probability P = P~( = 1 — P 2 £ and the two phases (£, rj) must be calculated numerically, 
by implementing the MSW Hamiltonian of the form given in Eq. (|A6|) . with number densities N e (x) and N n (x) 
evaluated along the neutrino trajectory in the Earth interior. From Eq. (|C2[) one obtains for the transition probabilities 
into electron neutrinos 
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Pel = \A n Zn+A 12 Z; 

Pe2 = \A\i Z\2 + A\2 Z22 
Pe3 = A 13 = Ue 3 , 



P- 



A 1 



JJ 2 

u e4 ' 



(C4) 
(C5) 
(C6) 
(C7) 



which, taking into account the expressions of the elements of the first row of the matrix A (An, A12, A13, A\± = 
ci 4 ci3,0, ci 4 si 3 ,si 4 ), become 

(C8) 
(C9) 
(CIO) 
(Cll) 

which constitute the generalization of the 3-flavor expressions used in the literature. It should be noted that in our 
parameterization it is Ay 2 = 0, and we can express the 4v transition probabilities in terms of the 2-flavor transition 
probability Pi", the knowledge of the phases (£,??) being unnecessary. For the transition probabilities into sterile 
neutrinos, one obtains the analogous expressions 



Pel 


_ 2 2 p2u 
— c 14 t 13- t el 


Pe2 


_ 2 2 p2v 


Pe3 


2 2 
— c 14 s 13 ! 


Pe4 


— s 14 ) 



P sl = \ A 41 Z 11 +A± 2 Z 2 i\ 



Ps2 
Ps3 
Ps4 



A41 Z\2 

u 



Aao Zo 



A 2 = 

A 43 — 

A 2 = 

A 44 — 



u 



s3 - 
2 

s4 ' 



(C12) 
(C13) 
(C14) 
(C15) 



which, in the limit of no matter effects (Z12 = R12), return the expressions of given in Eqs. (|1HI14[) . Now, differently 
form the case of the electron neutrinos, both elements A41 and A42 can be different from zero [see Eqs. (|15M18|) ]. and 
one cannot express the 4-flavor transition probabilities P s \ and P S 2 in terms of the 2-flavor transition probability 
P 2 e, needing the complete information on the (complex) 2-flavor transition amplitudes, having phases (£,fy). We 
close this appendix by stressing that the expressions given in Eqs. (|C4IC7[) and Eqs. (|C12IC15[) are valid for any 
parameterization of the mixing matrix of the form U = ARi2- The validity of Eqs. (|C8IC1 1| ) is instead restricted to 
the case in which the matrix A is of the form A = BR14R13. 



Appendix D: Potential sensitivity to the CP violating phases 



For definiteness, in the rest of the paper we have restricted ourselves to the case of vanishing CP violating phases. 
Here, we briefly comment on the potential sensitivity of solar neutrinos to them. In the 3-flavor framework one can 
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always eliminate the CP violating phase S appearing in the mixing matrix from all the observable quantities involved 
in the solar neutrino transitions [77I . l78j , since the two following conditions hold: I) All the relevant information on the 
flavor conversion is contained in the survival probability P ee of the electron neutrinos, as one cannot distinguish 
from v T at low energies and by unitarity it is P epi +P eT = 1 — P ee ] II) The phase S can be eliminated from the expression 
of P P p since it can be rotated away from the MSW dynamics due to the particular form of the potential [z3, [lH (see 
also [zISqI). 

The same conclusion is not true in the 4-flavor case since the two conditions above are no more valid, as we briefly 
show. Without loss of generality we can assign one of the three phases to the (2, 3) sector since, as we have seen in 
Sec. II and Appendix A, the associated mixing angle #23 can be eliminated from the description of the solar neutrino 
conversion. We can then attribute the remaining two phases to the (1,3) and (1,4) sectors, by defining the complex 
mixing matrix as 

U = i?23^24-R34-Rl4-Rl3-Rl2 , (Dl) 

where Rij is a 4 x 4 complex rotation in the plane, formed by replacing the real 2x2 submatrix in Eq. ([2]) with 
the complex one 

J.). (D2) 
\—Sij c tJ j 

with £ij = cosOjj, Sij = sin #jj e ll5ij . The two phases $13, $14 will appear in the expression of the transition probability 
P es in Eq. (|A23[) through the (now complex) elements U S i [see Eqs. (|1HI14[) ]. Furthermore, they will enter at the 
dynamical level, by affecting both P ee and P es , through the expressions of the mixing elements U™ in matter, which 
are determined by the diagonalization of the Hamiltonian 

frdyn_ v lnA (l 2 +r x \u\ 2 r x a*P\ , . 

where the two parameters a and /3 can be complex numbers. With our choice of the phases, /? = — S24 is still real, 
while a is complex and is obtained with the replacements (S13 — > §13, Su — > §14) in its expression in Eq. (|A7|) . 

We postpone to a future work the study of the complex extension of the treatment we have provided for the real 
case. Here we just limit ourselves to observe that the CP phases always appear in terms involving two small mixing 
angles, and therefore it is difficult to observe their effects in current solar neutrino experiments. Finally, we stress 
that for #24 = #34 = the CP phases disappear from the description of the solar neutrino transitions, thus rendering 
the numerical results presented in Sec. Ill independent of them. 
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